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In this talk I review the recent progress made in the calculations of quarkonia production 
in fixed target experiments. NRQCD organizes the calculations in a systematic expansion 
in ot s and v, the relative velocity between the heavy quarks. Within this formalism 
there are octet contributions which are not included in the color singlet model. These 
contributions depend upon unknown matrix elements of local operators which are fit 
to the data. Using these fits, there are several predictions which do indeed improve 
agreement with the data. However, the prediction for the polarization of the produced 
states as well as the ratio of the xi to X2 cross sections differ substantially from the 
data for the case of pion beams. Possible large corrections from higher twist effects are 
discussed as is the issue of the the proper choice of masses. 



1. Introduction 

The use of non-relativistic QCD (NRQCD) allows us to calculate both produc- 
tion and annihilation rates of heavy quark bound states in a systematic expansion 
in a s and v, the relative velocity of the heavy quarks!. Moreover, the inclusion of 
higher Fock states, which emerge naturally in the formalism, allows for a consistent 
factorization of long and short distant effects, thus validating the use perturbative 
QCD to calculate the Wilson coefficients. Furthermore, the long distance effects 
are now written in terms of well defined operators, which can be calculated on the 
lattice i, instead of potential model wave functions. 

In light of this progress, it is interesting to revisit iffl the issue of hadro- 
production in fixed target experiments. Previous calculations within the confines of 
the color singlet model B were found to be inconsistent with the datai. The predic- 
tion for the overall normalization of the cross section is too small as is the ratio of 
the production cross sections for \i an d Xi- The direct production rate for J ftp is 
too small, and the J/ip are predicted to be partially transversely polarized, which 
they are not. This disparity between theory and data is, at least for the first two 
observables discussed, crying out for a new production channel. NRQCD supplies 
just such a channel, namely the state in which the two heavy quarks are in a relative 
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octet configuration. This state will have a finite overlap after soft gluon emission 
acts as a color sink. The overlap will be suppressed by powers of v as dictated by 
the velocity scaling rules 13. 

Before proceeding to the results I would like to briefly discuss the levels of 
rigor which go into the various approximations in the calculations. First, there 
is no operator product expansion in these calculations. Thus, the factorization is 
performed via a diagrammatic analysis, as is done for Drell-Yan and other such 
processes. Factorization in such cases in known to be violated by higher twist 
effects which are suppressed by powers of the large invariant mass scale involved in 
the process. For the case discussed here, this scale would be quarkonium mass. In 
this regard, the proofs of factorization in the case of small px production is no less 
rigorous than at large pr- The only difference in the two cases is that the higher 
twist corrections at large pt are suppressed by l/p T as opposed to l/m? H . Thus, 
we expect larger errors to be incurred at small px- 

In showing factorization it is imperative that one sum over all relative color 
states of thequarkonium, otherwise, as was discussed in the case of the decay of P 
wave states there will be an infrared divergent Wilson coefficient which obviously 
destroys any hope of calculating in a model independent fashion. Once the higher 
Fock states are included, the factorization is restored and the final result of our 
calculation depends upon unknown non-perturbative matrix elements which are 
enumerated by the velocity scaling rules. Thus, the production cross section for the 
reactions 

A + B — > H + X, (1) 

can be written as 

1 

a H = Y1 J dx l dx 2 fi/A(Xl)fj/ B (X2) <j{ij -> H) , (2) 

i,j o 

^3^ H ) = E C QQln ] (°n)- (3) 
n 

Here the first sum extends over all partons in the colliding hadrons and fi/A etc. 
denote the corresponding distribution functions. The short-distance (x ~ 3> 
l/{m,Qv)) coefficients C^qw describe the production of a quark-antiquark pair in 
a state n and have expansions in a s (2mg). The parameters (O^) describe the 
subsequent hadronization of the QQ pair into a jet containing the quarkonium H 
and light hadrons. 

The velocity scaling rules are derived via the multipole expansion, which tells 
us that soft gluon couplings to a heavy quark bound state are suppressed by the 
ratio of the size of the bound state to the wavelength of the gluon. This expansion 
was first used within the confines of the strong interaction by Gottfried t2l a while 
back, and it seems to work quite well. Furthermore, present extractions of matrix 

"Their precise definition is given in Sect. VI of §. 
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elements seem to agree with their predicted scalings with v. As such, we will 
continue under the assumption that these scaling rules are valid. Furthermore, we 
will ignore effects due to the finite size of the target. Presumably, the soft gluons 
which will be exchanged between the target and the quark-antiquark pair during 
the hadronization process should lead to higher order effects in v 2 for the same 
reasons as stated above. We will assume this is true for now, and will keep this, 
perhaps dubious assumption, in the back of our minds when we confront the data. 

2. ip' production 

The case of ip' production is simplest to analyze since there are no states below 
open charm threshold which contribute to its indirect production rate. Many of the 
arguments discussed in this simple case will apply for the other states as well. 

The production of ip' in the singlet channel begins at 0((^) 3 w 3 ) due to charge 
conjugation, while the octet channel production is 0((^) 2 v 7 ). Given that numer- 
ically, a (4m;;) oc v 2 , it would seem that octet production should be of the same 
order as the singlet channel. However, the singlet cross section vanishes at thresh- 
old where there is small x enhancement due to the gluon distribution functions and, 
as such, the octet actually dominates the singlet. 

Before going on to quantitative issues however, we must address the issue of the 
proper choice of the hadron mass. It is clear that the short distance coefficients 
should be calculated using the quark mass, since binding effects are neglected by 
definition in this part of the calculation. However, at face value it seems that as 
far as the phase space boundary is concerned we should be using the hadron mass 
instead of twice the quark mass. This issue was vehemently debated during the 
workshop. As I emphasized then, the proper choice of mass in the phase space 
boundary is indeed the 2m c . This is the choice which is consistent with the v 2 ex- 
pansion, as the binding effects are always higher order in this expansion parameter. 
This is best illustrated in the case of heavy-light meson decay where the expansion 
parameter is Aqcd/toq. In this case, we may perform an operator product ex- 
pansion, with no question as to the proper choice of mass. At leading order in the 
OPE, we find that the phase space is dictated unambiguously by the quark mass 
as a consequence of unitarity. If we consider the lepton spectrum, then there is 
an explicit factor of 6(1 — 2Ei/itiq) in the differential rate. The fact that the true 
phase space boundary is determined by the meson mass is seen when one goes to 
higher order in the OPE where the expansion looks like 

oc 0(1 - 2E l /m Q ) + 6(1 - 2E l /m Q ) + 6'(1 - 2E l /m Q ) + . . . . (4) 

We see that the expansion breaks down near the partonic endpoint Ei = 2toq and 
is signaling the need for a resummation of the non-perturbative effects. Such a 
resummation leads to the construction of a structure function which has support all 
the way to the hadronic endpoint. In our case a resummation of higher order v 2 ef- 
fects will lead to shifting the space limits from being partonic to hadronic. Whether 
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or not the non-perturbative corrections are large depends upon the observable of 
interest. In the example above we may safely use the partonic mass if we are not 
interested in the endpoint region. With that said, let us confront the theory with 
the data. 

In previous analyses, performed using the color singlet model, it was found that 
the theoretical predictions needed a K factor of 25 a. However, this discrepancy is 
greatly reduced once we use twice the charm quark mass instead of the hadron mass. 
Indeed, given the uncertainty in the charmed quark mass, varying m c between 1.3- 
1.7 GeV, changes the total cross section by a factor of 8 at */s = 30 GeV. This 
large variation is a consequence of the steep rise in the gluon distribution at small 
x. Nonetheless, let us press on assuming that the expansion is well behaved and 
consider the consequences. In ref i, it was found that, using m,H = 2m c = 3 GeV 
the color singlet contribution fell a factor of 3 below the data. Including the color 
octet contribution leads to a fit of the data with the choice Ag = 5.2 • 10 -3 GeV . 
At large p T , the authors of found (Of { y S Q )) + (Of ( 3 P )) = 1.8 -lO" 2 GeV 3 . 

If we assume (Of ( 1 S'o)) — (£2f ( 3 Po)} /m 2 ., the fixed target value is a factor of four 
smaller than those found in til. This discrepancy should not concern us, nor should 
we consider this particular observable a good test of the color octet mechanism for 
the reasons discussed above. 



3. J/tp production 

Using the data from proton beam fixed target experiments we may again fit the 
data using the value 

Ag(J/ip) = 3.0- 1CT 2 . (5) 

As in the previous case this observable is very sensitive to the choice of the quark 
mass and, as such, the fact that the color singlet contribution is below the data is 
not strong evidence for the existence of the octet channel. However, in this case we 
may also look at the ratio for the direct to total cross section which is not sensitive 
to the quark mass. Indeed, we find that the pure singlet contribution gives a ratio 
0.21 whereas inclusion of the octet with the matrix element extracted above givesl 
0.63 which is in much better agreement with the experimentally extracted value of 
.62 ± 0.04 for a proton beam at i/s — 23.7 GeV t3. Note that this is not a trivial 
consequence of fitting the color octet matrix element since the indirect contribution 
is dominated by color singlet gluon fusion and the singlet matrix elements are fixed 
in terms of wavefunctions E3. 

Another very interesting observable is a Xl /o X2 which has been measured in 
proton as well as pion beam experiments. The singlet cross section for \i production 
is suppressed by a factor of a relative to \2 cross section, while the leading octet 
contribution is (0(a 2 v 3 ) but is suppressed because it is a quark initiated process. \2 
production on the other hand, is dominated by color singlet gluon fusion at 0(a 2 v 5 ). 



The E771 experiment measured a value 



(T X1 /<J X 2 = 0.34 ±0.16 for a proton beam 
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Fig. 1. Total (solid) and singlet only (dotted) tp' production cross section in pion-nucleon collisions 
(xp > only). The solid line is obtained with Ag(ip') = 5.2 ■ 10 — 3 GeV 3 determined using the 
data from proton beam experiments. 



at y/s = 38.7 GeV, and NRQCD predicts a value of O.OTu? However, the relativistic 
corrections can be substantially given that the leading singlet contribution scales like 
a 3 v 3 and numerically a(2m c )/ir oc v 3 . Indeed at 0(a 2 v 9 ) there will contributions 
coming from intermediate 1 So and 3 Pj octet as well as singlet states, as well a octet 
3 Dj states. Though these states are suppressed in v 2 given the large number of 
channels which contribute, the net contribution could be substantial. A naive use 
of the velocity scaling rules leads to ~ 0.3? 

For pion beams theory does not seem to do as well. The latest reported value 
for this ratio for pion beams is given by 0.57 ± 0.18, whereas the leading NRQCD 
prediction is given by 0.07 even including the next order corrections in v 2 , it is clear 
that the theory falls short. Furthermore, the overall normalization for the total J/vE' 
and tp' production cross sections, found using the fitted values of the octet matrix 
elements using the data from proton beam experiments, also falls short as is shown 
in figures 1 and 2? 

4. Polarization in fixed target experiments 

There are some interesting theoretical issues involving polarized cross sections 

6 In ref. E the ratio was weighted by the branching ratios to J/ifi. 

c The value quoted in the first reference in Ej was not weighted by branching fractions of \J m t° 
J/tp, contrary to what is stated in the text, m 

d data points in the plot for figure 2 in ref. Q are off set in y/s due to an error in the plotting 
routine. 
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Fig. 2. J/ip production cross sections in pion-nucleon collisions for xp > 0. Direct J/ip production 
in the CSM (dashed line) and after inclusion of color-octet processes (dotted line). The total cross 
section (solid line) includes radiative feed-down from the \ c j and ip' states. The solid line is 
obtained with A s (J/ip) = 3.0 • 10" 2 GeV 3 . 



within the NRQCD formalism. However, due ±o space limitations, they will not be 
discussed here, and I refer the reader to refs. Ii3tht3 for discussions. As we will see 
polarized production is a useful tool for investigating the octet mechanism. 

Polarization measurements have been performed for both ip and ip' @ pro- 
duction in pion scattering fixed target experiments. Both experiments observe an 
essentially flat angular distribution in the decay ip — > [ip = J/ip, ip')i 

cx 1 + a cos 2 9 , (6) 

dcoso 

where the angle 8 is defined as the angle between the three-momentum vector of the 
positively charged muon and the beam axis in the rest frame of the quarkonium. 
The observed values for a are 0.02 ± 0.14 for ip' , measured at y/s — 21.8 GeV in 
the region xf > 0.25 and 0.028 ± 0.004 for J/ip measured at y/s = 15.3 GeV in the 
region xp > 0. In the CSM, the J/^'s are predicted to be significantly transversely 
polarized □, in conflict with experiment. 

The polarization yield of color octet processes can be calculated along the lines 
of the previous subsection. We first concentrate on ip' production and define £ as 
the fraction of longitudinally polarized ip' . It is related to a by 
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For the different intermediate quark-antiquark states we find the following ratios of 
longitudinal to transverse quarkonia: 





1 : 3.35 


£ = 0.23 


lc(8) 
^0 


1 : 2 


t = 1/3 


3p(8) 


1 : 6 


£ = 1/7 


3 S (8) 


: 1 


e = o 



(8) 



where the number for the singlet process (first line) has been taken from Ef Let us 
add the following remarks: 

(i) The 3 S' 1 8 ^-subprocess yields pure transverse polarization. Its contribution to 
the total polarization is not large, because gluon-gluon fusion dominates the total 
rate. 

(ii) For the 3 pj 8 ' ) -subprocess J is not specified, because interference between 
intermediate states with different J could occur as discussed in the previous sub- 
section. As it turns out, interference does in fact not occur at leading order in a s , 
because the only non-vanishing short-distance amplitudes in the JJ Z basis are 00, 
22 and 2(— 2), which do not interfere. 

(iii) The ^subpr ocess yields unpolarized quarkonia. This follows from the 
fact that the NRQCD matrix element is 



(0tfT A ^a$\$ ^T A X \Q) = 1 (Of'CSo)) , (9) 

independent of the helicity state A. At this point, we differ from 0, who assume 
that this channel results in pure transverse polarization, because the gluon in the 
chromomagnetic dipole transition 1 Sq S " > — ► 3 s[ 8 ^ + g is assumed to be transverse. 
However, one should keep in mind that the soft gluon is off-shell and interacts 
with other partons with unit probability prior to hadronization. The NRQCD 
formalism applies only to inclusive quarkonium production. Eq. (^) then follows 
from rotational invariance. 

(iv) Since the 3 Pj and 1 Sq -subprocesses give different longitudinal polariza- 
tion fractions, the ip' polarization depends on a combination of the matrix elements 
(Of^So)) and (Of ( 3 P )) which is different from A 8 (V')- 

To obtain the total polarization the various subprocesses have to be weighted 
by their partial cross sections. We define 



58{H) ~ Vffj (10) 



and obtain 



e This number is xp -dependent and we have approximated it by a constant at low Xp, where the 
bulk data is obtained from. The polarization fractions for the octet 2 — + 2 parton processes are 
^-independent. 



7 



0.23 



„ f3o( 1 )\ 



i* 8 (^) + I(l-tf 8 (V,')) 



= 0.16 + 0.11^(^0, 



(11) 



where the last line holds at y/s — 21.8 GeV (The energy dependence is mild and 
the above formula can be used with little error even at y/s = 40GcV). Since < 
5 8 (H) < 1, we have 0.16 < £ < 0.27 and therefore 



0.15<a<0.44. (12) 

In quoting this range we do not attempt an estimate of 5s(ip'). Note that taking 
the Tevatron and fixed target extractions of certain (and different) combinations 
of (Og ( 1 5o)) and (Of ( 3 i"b)) seriously (see Sect. 5.1), a large value of 8g(ip') and 
therefore low a would be favored. Within large errors, such a scenario could be 
considered consistent with the measurement quoted earlier. From a theoretical 
point of view, however, the numerical violation of velocity counting rules implied 
by this scenario would be rather disturbing. 

In contrast, the more accurate measurement of polarization for J ftp leads to a 
clear discrepancy with theory. In this case, we have to incorporate the polarization 
inherited from decays of the higher charmonium states XcJ and ip'. This task 
is simplified by observing that the contribution from XcQ and Xd feed-down is 
(theoretically) small as is the octet contribution to the Xc2 production cross section. 
On the other hand, the gluon-gluon fusion process produces Xci states only in a 
helicity ±2 level, so that the J/ip in the subsequent radiative decay is completely 
transversely polarized. Weighting all subprocesses by their partial cross section and 
neglecting the small ip' feed-down, we arrive at 



£ = 0.10 + 0. 116 8 (J/ip) (13) 

at y/s = 15.3 GeV, again with mild energy dependence. This translates into sizeable 
transverse polarization 



0.31<a<0.63. (14) 

The discrepancy with data could be ameliorated if the observed number of Xci from 
feed-down were used instead of the theoretical value. However, we do not know the 
polarization yield of whatever mechanism is responsible for copious Xci production. 

Thus, color octet mechanisms do not help to solve the polarization problem 
and one has to invoke a significant higher- twist contribution as discussed in B. 
To our knowledge, no specific mechanism has yet been proposed that would yield 
predominantly longitudinally polarized ip' and J ftp in the low xp region which 
dominates the total production cross section. One might speculate that both the 
low XcilXci ratio and the large transverse polarization follow from the assumption 



8 



of transverse gluons in the gluon-gluon fusion process, as inherent to the leading- 
twist approximation. If gluons in the proton and pion have large intrinsic transverse 
momentum, as suggested by the pt-spectrum in open charm production, one would 
be naturally led to higher-twist effects that obviate the helicity constraint on on- 
shell gluons. 

5. What Needs to Be Done 

The uncertainties in the theoretical prediction at fixed target energies are sub- 
stantial and preclude a straightforward test of universality of color octet matrix 
elements by comparison with quarkonium production at large transverse momen- 
tum. Small-rz;, as well as kinematic effects, could bias the extraction of these matrix 
elements in different directions at fixed target and collider energies!. The large 
uncertainties involved, especially due to the charm quark mass, could hardly be 
eliminated by a laborious calculation of a s -corrcctions to the production processes 
considered here. To more firmly establish existence of the octet mechanism there 
are several experimental measurements which need to be performed. Data on po- 
larization is presently only available for charmonium production in pion-induced 
collisions. A polarization measurement for a proton beam would be very interesting 
given that we seem to have a better handle on the theory in this case, as is demon- 
strated by the observed value of the X1/X2 ratio. A measurement of polarization at 
large transverse momentum or for bottomonium is of crucial importance, because 
higher twist effects should be suppressed. Furthermore, a measurement of direct 
and indirect production fractions in the bottom system would provide further con- 
firmation of the color octet picture and constrain the color octet matrix elements 
for bottomonium. 

From a theoretical standpoint there are still several issues that warrant further 
investigation. To begin with, factorization in hadro-production of quarkonia is 
presently just a working hypothesis. This is true at large pr as well as at small 
Pt- Of course while formally this puts both calculation on the same theoretical 
footing; practically, there is still an important difference between the two cases. 
Namely, the size of the higher twist effects will be suppressed by 1 jp\ at large px , 
as opposed to l/(4m^), for our fixed target cross section. It is furthermore possible, 
and this again applies to the case of large pt as well, that the higher twist effects 
could be enhanced by powers of \/v. Indeed, understanding higher twist effects in 
quarkonium production is complicated by the presence of the scales mv and mv 2 . 
Needless to say, there is still much work to be done on this subject before we can 
get a good handle on the errors due to higher twist effects. 
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